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ABSTRACT 

The present work is concerned with the study of thermal stress analysis of a two-dimensional rectangular plate under the 
hyperbolic heat conduction model with a heat source inside the body under the constant surface temperature. The 
internal heat source varying periodically with time is considered. The thermoelastic displacement, stresses and 
temperature are determined within the context of non-Fourier heat conduction theory. The solution for temperature field 
and stress field are obtained with the help of differential transform method and later used to find displacement. The 
effect of periodically varying heat source on thermal stresses is studied. The special case is considered and results are 
illustrated numerically and graphically for a copper plate. 
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INTRODUCTION 

The study of the hyperbolic heat conduction equation in the field of thermoelasticity plays a vital role in many 
branches of science and technology. Thermoelasticity is the study of deformation occurs in the body due to thermal 
and mechanical load applied on it. As the advancement of technology, in many experimental results, it is observed 
that, heat is conducted in a material with a finite speed which contradicts the fact of infinite heat of propagation 
which is explained well by Fourier law of heat conduction or parabolic heat conduction. The heat is transported in 
hyperbolic form when the temperature gradient is very high, extremely large heat fluxes and time duration is very 
short. Cattaneo [1] and Vernotte [2] proposed a modification in the classical theory of heat conduction by bringing a 
new term called relaxation time for heat flux in solids, which is given by the relation 

f?+^g = -ivr (i) 

where, tv- is called relaxation time, tj is the heat flux and Jl is the thermal conductivity of the material. 

Most of the research works on the thermal response of rectangular plate in the parabolic heat conduction 
domain are available in the literature. V. R. Manthane et.al [3] obtained the solution of the transient thermoelastic 
problem of non-homogenous rectangular plate using integral transform technique. N. Sumi [4] investigated the 
thermal stress in a finite rectangular plate by the complex method. T. C. Lim [5] investigated the thermal stress in a 
thin auxetic plate and provides a different choice to decrease the thermal stress. Chen [6] has studied the thermal 
stress in a rectangular plate subjected to non- uniform heat transfer in terms of power series using Lanczos- 
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Chebyshev and discrete least square methods. Tanigawa and Komatsubara [7] have discussed the thermal stress and its 
intensity factor under compressive stress field in a rectangular plate. From the last few decades, the researcher works on 
non-Fourier heat conduction. As per El Dhaba and Abou-Dina [8] analysis, the thermal stress in a long cylinder of 
rectangular cross section occurs due to variable heat sources and variable pressure at its boundary in finite Fourier 
transforms domain. Chen and Hu [9] suggested that HHC model are more efficient to design structure against fracture 
under thermal loading. Wang and Li [10] studied the HHCE for the fracture of a finite piezoelectric layer using Laplace 
and dual integral transform and also investigated thermal shock resistance in a plate under HHC model. Ahmed [11] 
calculated the thermal stresses numerically for a nano composite slab. M. R. Talaee and V. Sarafrazi [12] presented an 
analytical solution for non Fourier heat conduction in a solid of rectangular shape using Eigen function expansion and 
shows the use of laser absorption in tissues. H. R. Al-Duhaim et.al [13] determined the temperature distribution and 
thermal stresses for the laser short pulse heating of solid surface considering HHC model. 

Yang [14] has applied a sequential method to find the solution of direct and inverse HHCE using finite difference 
and the Newton Raphson method. N. S. Al-Huniti and M. A.Al-N.M.Naji [15] worked on the thermoelastic nature of rod 
under the effect of moving heat source considering HHC model. L. Ping et al.[16] worked on the higher-order 
thermoelastic response to determine the thermal bending of the laminated composite plate. M. H. Hsu [17] presented the 
solution of the hyperbolic heat conduction equation and also investigated the impact of relaxation time on temperature 
distribution. I.V.kudinov and V.A.kudinov [18] studied the dynamic stresses in an infinite plate associated with HHCE. 
B.A.Hamid [19] worked on an analytical solution of non-Fourier heat conduction equation with periodic thermal 
oscillation using integral transform technique. In past literature, many authors used DTM (Differential Transform Method) 
to solve different initial and boundary value problems and showed the result obtained is much closer to the exact solution. 
Afshin Babaei et.al [20] discussed the inverse problem of heat conduction using reduced DTM. Shaher momani and Zaid 
odibat [21] discussed the solution of non-linear fractional partial differential equation applying DTM. Jafari, Sadeghi and 
Biswas [22] have studied an analytical solution of multi-dimensional partial differential equation by DTM. This method, 
due to its simplicity and effectiveness are used by many researchers. 

The present work is concerned with the study of thermal stress analysis of a two-dimensional rectangular plate 
under the hyperbolic heat conduction model with a periodically varying heat source inside the body under the constant 
surface temperature. The thermoelastic displacement, stresses and temperature are determined within the context of the 
generalized theory of thermoelasticity. The solution for temperature field and stress field are obtained with the help of 
differential transform method. The effect of periodically varying heat source on induced thermal stresses in a plate is 
illustrated numerically and graphically for a copper plate. 

MATHEMATICAL MODELLING LOR TEMPERATURE FIELD 


Consider a plate of rectangular shape occupying the space D: [(sr.. y, t ) : x E [0, a], y E [0, i], t E [O. T]]. The unsteady state 
temperature distribution under the hyperbolic heat conduction in the rectangular plate with heat generation inside the body 
is given by 


<tT a*r_ /a- a- \ p 
“77 + T., — - k + 77F + 77 


( 2 ) 


where, k = — ,p,c are the thermal diffusivity, mass density and specific heat capacity respectively and Q is the internal 
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heat source. 

With initial and boundary conditions are given by 


tUi v, g) = r B! 

(3) 


(4) 

7(0., y, i) = 7, 

(5) 

t£x, a, t) = 0 

(6) 

Tts, y, i) = 0 

(7) 

Tis, b, t) = 0 

(8) 


Solution for Temperature Field 


Considering T 0 = and T, = ‘V sin (~[r) ain (~) 

And solving of the problem defined from eq.(2) to eq.(8) is done semi analytically using differential transform 
method given in reference [16]. Thus, after applying DTM on eq.(2) to eq.(8), one will get 


0 + 2) = r ( a + £,| 8+ ,-, [W + lXf + ZMf +■ 2, ri , 03 + irfrj + 1) (fj + 2)7(<+2.9)- (( + >}. 0 4-1) + 

SffiSai 

P€ \ 


T(f. f ,.o) =^ =n=i3 ,^ m ' nri ":;: L ' r =o 




(9) 




~rr fr, rV'« TPS l&T® m l n H ^‘ l -ll< ■'Z 7 n Q, 'i n 

= s n=Mi , D=liil ^- — -, nf..o.0) - 0 


( 10 ) 


Now, using these transformed boundary and initial conditions in eq.(9), one can obtain 7(< , T), 9) for different 
values of rj, 0 and which is given as: 




m.=c=i,3.5,, 


^ 4- 4- 


^rCc^-om 


F=0/2 e ,.!T 4| t * 


(ID 


where, f, tj are odd and 9is even 


7C£,jj,0) = X£ 


: 12=1,3, i,„ KBJr Z 


^rCCfj.ojI 


B(-lJP“ 1 kPTqP 


F =^r i 0 - 


(urir 4- n 2 Tt :: 3 p 4- 




( 12 ) 


where, ij, 9 are odd 

Now, we considered here the heat source in rectangular co-ordinates Q(x, y, t) which varies periodically in the 
following form: 
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f?U, y, f) 


Sfjc)GCy} 
} ah 



0 < t < U 


= 0 t > T 

Where^j., the strength of the heat source and 5 is the Dirac’s delta function. 
On applying differential transform on equation (13), we get 


(13) 






aa 0 


r sm l77 


(14) 


By applying inverse differential transform and using equation (11), (12) and (14), one can obtain series solution 
for T(x,y, t) which is as follows: 


Tfc y, i) = rf =0 2^ =0 2e =0 n<^, 9)*Vi e (15) 

Using 0) for each values of £, f), 0 from eq. (11), (12) and using (14) in eq. (15), we can get the solution 

for temperature function as: 


TU.y ti JfflLdn(S=.W'HEJ - -, J*2L. b AS)(!S_ rfa S 

"-‘■••“own 2 - a / ■ c t rrftW Vr r ■■ titi - 

r= iS7 a . /war ’i , 2 . . ,» * - , z I’l 

**£=»=IAS.»«■(. « j 5 - u l. s J taW + rfwVilfcs 9. + | 


MATHEMATICAL MODELLING LOR THERMAL STRESSES 


06) 


Considering thermal stress function % in rectangular coordinates system defined in reference [23], as given by 


X = Xc +■ Xp ■ 

Where % c is a complementary function and is a particular integral. 
And find are satisfied by the equations 



And component of stress function are given by 


n -^ X 

09) 

a - 

(20) 

~ 3x3y 

(21) 
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And the boundary conditions 

a** = Oi Gjfj, =0 on s = ft o (22) 

Solution for Stress Function 

The complementary function % c can be obtained by solving eq.(17). 

Let + ^ + (23) 

Solution of particular integral jf_is obtained by applying differential transform on eq. (18). The transformed form 
of eq.(18) is given by 

(jh + 2 Jlgw -t- liF(pi + Z,v,w) + fv +■ 2)(i' + + 2,os) - S(u - f. v - if .ns - 0) (24) 


where, V (fi, T, a) is the transformed function of % 0 . 

For each fi r v r OJ> the values V(jl,v, of) by recursive method can be evaluated as: 


ET=^=i^x§SrC(t-l> r <* - 2 - 3*i**rTH - 2 - 2 t,t) + 2r.0>) + - 


2,t}, 0-1) 


Now applying the definition of inverse differential transform, one can obtain as: 

Xp = -aE Zf = ;Z^ =1 Zi =1 LE r ^^(-l3 r (< - 2 - Zr)!^^rCf - 2 - 2 r,ij + 2r,0) J.v f y r? i a + {”- 

S inf)x= 


(25) 


(26) 


Using eq. (23) and eq. (25), we will get expression for stress function^, which is given by 
* = ^+7 xy~ + fy* - «E IfU ^ E£ =1 L - 2 - 2 r) ! - 2 - 2 r., + 2 r. 9 )) + 

cE Q a fr Y firt . izt \ 2 

l7- sm 7J r (27) 

And now using this value of stress function in eq. (19) "*(21), we will get stress component as: 


= cjs* + W ~ aE I^a^=LSe=i.^S r i 0 ^C-l3 r (< - 2 - 2r)E^^r(( - 2 - 2r,r\ + 2 r, 


(28) 


e,y-rfI^I^SSU^s£(c-iy<l-I-ar)l^l1l-2-2r., + ar.e))»*^Vl» + aAg)*g- 

sin 7) 


(29) 
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a xy = -|« 1 *+ c zy- aE Sf=i^ , =i5*=i^E r i*(C-l5 , '<?-Z-2r)!^rr(f -Z-Zr.jj + 2r,e)j*f V-V + 

E£-^(l) S {E- aSn E)J 

13 TqP'C Vjt/ \t 7 / I 

Now applying the boundary condition defined in eq. (22), we will get 


(30) 


Cn — Cg — U 


c L = uE 2f = , 2f =1 iZ r z _ B ((-l3 r « - 2 — 2r3 ! (1 + 2r)l Kf - 2 - 2 r, 1 + 2r, 0 ])a?-- L j° - {7 - sin 7 ) aNow substituting 

these values of c i ,c Ii ,c g in eq. (28)-^(30), one will get the expression for stress components as: 


= -*ESf= s ^ =1 I^=4i: r i (f ((-ir« - 2-2-2-2r. T 42r.e>J*fy 


tf -2 4.0 


(31) 


O r Cf 2 2^1 <1 | 2r)!7« 2 2r. 1 I 2r, 9»af _1 4® {* 

sir. 7 I ay - aETf. s T.= Jg. L ^ ^ {c-lFtf - 2 - 2r) !^^T({ -2-2r.r+2 r 0 l) + 

- Z - Zr) I & 4 2r)I it? - 2 - 2r, 1 4 2r, - 7 ^( 7 )'{ 7 - 

+ aFSr =: S- =1 Se =1 ^S r f ff (<-l) r « - 2 - Zr)I^7« - 2 - 2r., + 2r.e)p-^ L t e - 

££_£L(l) : (IL_ sill £I). r 
13 r^pf a it/ Vi r > 

MATHEMATICAL MODELLING LOR DISPLACEMENT LUNCTION 

The displacement function defined in reference [23] for plane stress in a rectangular coordinate is given by 


U - —^ 

* seL ix i+v5v. 


? 2 fi I is L+ 1 ‘ 8 *. 


(32) 


(33) 


(34) 

(35) 


where, G and v are shear modulus of elasticity and poison’s ratio, and -p satisfy the following equation as 


tr m + dyy + sET - 77 ; and f ^-7 4 ^ 7 ) ^£r = 0 


3x3y Sxiy By 2 ! 

Solution for Displacement Function 


(36) 


&\$r _ &\jr 

To calculate displacement function given by eq. (34) & (35), firstly, one must calculate — & — which can be determined 
by taking integration on eq.(36) with respect to x & y respectively, therefore 
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dip 

Hx 


“fi ET= S 2 «= 1 ^ 2 £(C- D r Q: - 2 - Zr> ■ ft + -z - 2r, 1 + 2 r, 0 'i)cf~ L t & £ + (~ - 

■in")(* 7 - y) + (Ef=i^=i£SUItc.+ E£=, B>^fi e ) + 0 (x) 


(37) 


8 # . 
Sjr = 


-2-2r>!^r(t_2 —2.,-2r + + 

off Sf =: S9 = i^;S^((-l) r H' -2 - 2r)!(l + 2r)!r{< -2 - 2r. 11- 2r,S>)#-V*' - 7277(7^ (*- 
sin7} (c,x}’ - i) + aE('Zf-L^-L S9-3 Ti i.JI.0)77-1 0 + S^-1 Ssr(Q , 3.0)xy r ‘ t e } + 0(y} 


(38) 


Using the above values in eq. (34) & (35) and neglecting the rigid deformation, one can get the displacement 
function as: 


“r = 77[-2T. S Sf.xS|.^((-D r Cf-2 -2r)!^7tf-2-2r,, + 2r.9}) [(l + + 

_jl_] x S-z y «-H 6 + Sf = , If=i^ - 2 - 2 r) ! (1 + 2 r) lT(f- 2 - 2 r,l + 2 r, g - 

-^-1 _ StL (if (E_ a ; n E) f^_ -£-l! f -U - 

1+V 2 j T q pc Vir/ \ 7 T / L 2 1+V 2 1+vJ 

W - - if [Er= K ^^4£( C -l>^ - 2 - «|“H5 -2-lr., + 200 )} [(^) (^ + ^) - 

_ , _ 7 t .J :{ 1 + 7-3ETO: - 2 - 2 o 1 + [r-^] + 

-ffls -C r 'l = (—- stn—1 Co>r -1 ) [*-tH - (Sf= B) ^f S^ =; rCP.’?. e>^4)l 

rqprMr* 7 > ' L 1+vJ °— 7 ‘ fc+l l+v {| A — 1 I l+v/J 


Table 1: Properties of Copper 


a 

n x io-*f vr 


0 ,02s 

V 

0,33 

G 

4,5 X 10 Lf? 

E 

1.19 X 10 u N/m" 

k 

1.12S3 X 10 4 m 2 /s 






NUMERICAL RESULTS AND GRAPHICAL DISCUSSIONS 


(39) 


(40) 


To obtain the numerical results, the copper plate is used for which material parameter are taken from the reference [24] 
which are given in table 1. The stress variation along spatial direction is determined and plotted. All the values are 
calculated taking a = 1cm, b = 1cm, t = ls,T B = 0 and T, = 10 ’C 
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Figure 1: Stress Component XX Along X Direction. 



Figure 2: Stress Component YY Along Y Direction. 



Figure 3: Stress Component XX Along Y Direction. 
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Figure 4: Stress Component YY Along X Direction 



Figure 5: Stress Component XY Along X Direction. 



Figure 6: Stress Component XY Along Y Direction. 
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Figure 1 shows the variation of thermal stress component a xx , along x-direction. it seems that along x-axis 
compressive and tensile stress occurs due to the constant surface temperature. It is also observed from the figure that 
compressive stress occurs only between 0.7< x < 0,9 and in remaining part tensile stress is observed. 

Figure 2 shows the variation of thermal stress component o vv , occurred in mid-point of the plate along y- 
direction. It seems that, initially no stress occurs in the plate and it decreases gradually as move along the outer edge of the 
plate as an effect of periodically varying internal heat source. It is also shown from the graph, as the time increases, the 
strength of stress increases in both the cases. 

Figure 3 shows the variation of thermal stress component a yy , along the vertical direction. It is observed from the 
figure that the values of stress component oscillate from compressive to tensile throughout along vertical direction. It is 
also clear from the figure that the nature of the stress component a xx is periodical. 

Figure 4 shows the variation of thermal stress component a vv , along horizontal direction. From the figure, it 
observed that initially there is no stress along x-direction but after the middle of the plate stress appears in a plate and 
continuously increases up to the end of the plate. 

Figure 5 shows the variation of thermal stress component , along the horizontal direction. From the figure, it 
is observed that varies linearly and shows maximum value at the boundary of the plate. 

Figure 6 shows the variation of thermal stress component a ^,, along the vertical direction. It is observed from the 
figure that the stress is constant throughout the plate along y-direction. 

CONCLUSIONS 

In the present work, thermal stress analysis of rectangular plate under the hyperbolic heat conduction model has been 
investigated. The response of the heat conduction equation under periodically varying heat source is also studied. The 
distribution of temperature, stress and displacement are obtained under the differential transform domain. From the 
numerical result and discussion, it is observed that 

• The stress component a xx is compressive and tensile in nature due to constant surface temperature and as an effect 

of relaxation time. 

• The stress component a vv shows opposite nature along spatial direction due to the periodically varying internal 

heat source 

• The stress component a xx shows different nature along x and y direction. 

• The strength of all the stress component increases as increase in time. 

From the present analysis, it is also observed that the periodically varying heat source affect the nature of the 
thermal stresses. And, also the internal heat source influences the temperature distribution under the hyperbolic heat 
conduction model. These results may be useful in the study of thermal characteristics of the various bodies in design and 
engineering applications. 
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